Numerical solution of the t-J model with random exchange couplings 

in (i = oo dimensions 
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To explore the nature of the metallic state near the transition to a Mott insulator we investigate 
the t-J model with random exchange interaction in d = cxd dimensions. A numerically exact solution 
is obtained by an extension of the continuous-time quantum Monte Carlo (CT-QMC) method to the 
case of a vector bosonic field coupled to a local spin. We show that the paramagnetic solution near 
the Mott insulator describes an incoherent metal with a residual moment, and that single-particle 
excitations produce an additional band, which is separated from the Mott-Hubbard band. 
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The extraordinary properties of high-Tj, cuprates are 
closely related to those of doped Mott insulators P, Q ■ 
In both systems strong electronic correlations play a key 
role. To understand their influence, fundamental elec- 
tronic correlation models such as the Hubbard and the 
t-J model have been studied intensively 0. In spite 
of their apparent simplicity these quantum mechanical 
many-particle models can only be solved approximately 
in dimensions d ~ 2,3. Thus the full range of physical 
phenomena described by the Hubbard or the t-J model is 
not yet understood, implying that their investigation still 
leads to unexpected, and often peculiar results. For ex- 
ample, it has been pointed that the Fermi-surface volume 
of the two-dimensional t-J model is inconsistent with the 
Luttinger-Ward theorem Indeed, by applying the 

Schwinger method to the t-J model Shastry [3] recently 
found a class of solutions which show precisely such a be- 
havior. In this situation it is desirable to obtain reliable 
conclusions about those correlation models at least in 
certain solvable non-trivial limits. In the case of doped 
Mott insulators "non-trivial" means that characteristic 
features such as a strong, but screened Coulomb repul- 
sion and the presence of local spin fluctuations in real 
space are retained. Namely, the screened Coulomb inter- 
action is responsible for the Mott metal-insulator transi- 
tion (MIT), and local spin fluctuations affect the quasi- 
particles by making the self-energy frequency dependent. 

The Mott MIT can be described by the dynamical 
mean-field theory (DMFT) 0,111 . The DMFT provides 
an exact, non-trivial solution of electronic lattice models 
with a local interaction such as the Hubbard model, in 
d = oo. It can be derived by mapping the original quan- 
tum lattice problem onto an effective quantum impurity 
coupled self-consistently to a dynamical fermionic mean- 
fleld ("bath") 0. However, due to the local nature of 
the DMFT inter-siic interactions are reduced to a static 
mean field. This implies that in the strong-coupling limit 
of the Hubbard model, which corresponds to the Heisen- 
berg model in the case of half-filling and to the t-J model 
in the doped case, non-local spin fluctuations are miss- 
ing. At the same time it is known from the investigation 



of spin models in the context of spin-glass problems, that 
when the spin coupling Jij is taken to be random non- 
local spin fluctuations survive even in d = oo, while the 
static mean fleld averages out [l0l - [T6| . In this case the 
self-consistency equations correspond to those of an ef- 
fective impurity problem coupled to a bosonic bath. In 
particular, the random-coupling Heisenberg model has a 
very remarkable property: as shown by Sachdev and Ye 
[l3| for SU(M) spins with M ~ oo its dynamical mag- 
netic susceptibility displays marginal Fermi liquid behav- 
ior, which was proposed in the phenomenological the- 
ory for the cuprate superconductors [T^]. This suggests 
a relation of the random-coupling Heisenberg model to 
the paramagnetic state of the cuprates. Indeed, spin- 
glass behavior was observed in La2_a;Sra;Cu04 in the low- 
doping regime a; = 0.04 [Tsj . 

The effect of doping on the random-coupling Heisen- 
berg model was studied in detail by ParcoUct and 
Georges [ll] in terms of the t-J model with random cou- 
plings Jij in the limit M ~ oo. This model also has a 
non-trivial d = oo limit and describes a doped Mott in- 
sulator without antiferromagnetic spin fluctuations. In 
particular, these authors calculated the coherence scale 
of quasiparticlcs and discussed the properties of an inco- 
herent metallic state near half-fllling. 

In this Letter we present a numerically exact solution 
of the random coupling t-J model for the realistic number 
of spin components, M = 2. This is made possible by an 
extension of the continuous-time quantum Monte-Carlo 
(CT-QMC) method to a local spin coupled to a bosonic 
fleld. We compute the quasi-particle energy scale as a 
function of doping and determine the spectrum of the 
incoherent metal near the Mott insulating state. 

The random coupling t-J model in d = oo. — To obtain 
a non-trivial d = oo limit (in the following we use the 
coordination number Z = oo instead) where non-local 
spin fluctuations are retained, the coupling constants Jij 
should be scaled as Jij = J*jl\fZ, with J*^ = const. [l9l - 
[2l| . in coinplete analogy to the scaling of the hopping 
amplitude However, since the static molecular fleld 

is proportional to Z J, the above scaling leads to a diver- 
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gence of the transition temperature for magnetic long- 
range order at half-filling, i.e., the paramagnetic state 
is unstable against an infinitesimally small external field 
for all finite temperatures. This problem does not occur 
when the coupling constants Jij are random variables, 
since then ^ Jij = but Ylj Jfj 7^ Oi implying that 
the static molecular field due to the surrounding sites 
averages to zero Doping of the random coupling 

Hcisenberg model then leads to the random coupling t-J 
model UM 



H = - 



'3 ' 



(1) 



where Ci^ = Cia{l - rii^^), Si = (l/2)X]cr<T' 
Here the summation is taken over nearest-neighbor sites, 
and we consider the Bethe lattice with infinite connec- 
tivity {Z = oo) p^ . The exchange couplings Jij are 
randomly distributed according to the probability distri- 
bution 



P{J,j) (xexp{-jfJ2J^ 



(2) 



Neglecting spin-glass order, the model defined by ([T]) 
with Z = oo reduces to an effective impurity model cor- 
responding to the action [l^ 



imp 



J drdr'l - /t(r)(?-i(r - r')/.(r') 



dTUnf^{T)nf^{i 



(3) 



where the repulsion U is taken to be infinite to exclude 
double occupation. The local electron is denoted by the 
Grassmann number fa, and we introduced nfa = flfa 
and Sf = (1/2) Y.aa' fl<^<y<y' f<y' ■ The local propagator Q 
and the local time-dependent exchange interaction are 
determined by the self-consistency conditions [l^ 



J{iVr. 



+ H- t^Gi,„p(ia;„), 



(4) 
(5) 



where Gimp is the single-particle Green function, and 
Ximp is the spin susceptibility, both evaluated in the ef- 
fective impurity model. Furthermore, lo„ and Vn are 
fermionic and bosonic Matsubara frequencies, respec- 
tively. The self-consistency condition corresponds to 
that of the extended DMFT j2SJp, M, HI with semi- 
circular density of couplings j26l l27j. 

The first and second term in 5'inip can be repre- 
sented by fermionic and vector bosonic baths, respec- 
tively, which we express by the operators a^a and bq. 
The corresponding Hamiltonian is written as 

ka q 

+ vY,{flaa + alfa)+gSr{b + b^), (6) 



where = N-'^/'^J^k^ka and b = N-^/'^J2qbq, with 
N as the number of sites. The quantities entering in ifimp 
are connected with those in S'imp by the relations 



k 



(7) 
(8) 



with Q^^{iijJn) = ii^n + M ^ A(ia;„). 

Spin-Boson Coupling in CT-QMC. — We solve the 
effective impurity model ^ using the hybridization- 
expansion solver of the CT-QMC [231 • An algorithm for 
the inclusion of a bosonic bath into the CT-QMC was 
formulated by Werner and Millis in the case of impu- 
rity models with electron- phonon coupling [2^ [sO] ■ They 
used the so-called Lang-Firsov transformation to elimi- 
nate the coupling term and thereby arrived at an effi- 
cient treatment of phonons. An exchange coupling can- 
not be eliminated in this way, since the spin operator Sf 
has three components which do not commute. For that 
reason wc treat the spin-boson coupling by a stochas- 
tic method. Namely, wc perform a double expansion in 
terms of the hybridization and the spin-boson coupling, 
and evaluate the series by Monte Carlo sampling. 

Results. — We now present the result of the model ([T]) 
as a function of the dimensionless coupling strength J/i, 
the particle density n, and the temperature T . In the 
following, energies are measured in units oiW = 2t= 1, 
where W is the half-width of the density of states of the 
non-interacting model. 

The quasiparticle renormalization factor z = [\ — 
ImS](za;o)/'^o] ~^ is shown in Fig. [ija) as a function of 
density n. While for J = 0, z tends to zero for ?i = 1, 
for J /t = 0.5 it approaches zero at n ~ 0.66 = ??o with 
increasing slope. This behavior is contrast to the result 
in M = oo, where z vanishes at n = 1 regardless of the 
value of J [l^. In the realistic case M = 2 the system 
hence appears to be a non- Fermi liquid for n > uq. 

To identify the type of electronic state which is stable 
for n > no we calculate the residual moment defined by 
x(t = /?/2), which is equal to Tx in the limit T — J- 0, but 
which approaches the T ~ Q value faster than Tx- The 
residual moment per electron is plotted in Fig.[TJb). The 
value at n = 1 corresponds to 6*^/3, which is consistent 
with previous QMC calculations [ij]. We see that x(/3/2) 



increases with decreasing T for n 



while it decreases 



for n < uq. From this data we conclude that the phase 
has a residual moment for n > hq. 

Next, we compute the momentum- resolved single- 
particle excitation spectrum A(e,uj), which is defined 
in terms of the local self-energy = Q~^{iu!n) — 

G;TJ^(ia;„) by 
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fi — e — E(a;+) 



(9) 



Here e = e^, with e[— W, W], is the dispersion in d = oo 
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FIG. 2: Intensity plots of the single-particle excitation spec- 
trum A{e,ijj), and the integrated intensity A{ijj) for J/t=0.5, 
n=0.6, 0.8 and T=0.0025, 0.025. The black solid curve is the 
analytical continuation from the spin- averaged G(iaj„), while 
red dashed lines correspond to those for individual spin com- 
ponents. 



FIG. 1: Three physical quantities characterizing the param- 
agnetic state are plotted as a function of n: (a) the renormal- 
ization factor z, (b) the residual moment x(/3/2) per electron 
normalized by = 12, and (c) the effective chemical po- 

tential /ieff = M ~ ReE(O). The dashed curves in (c) show the 
chemical potential of the paramagnetic and polarized state, 
fipM and /iFM, for a noninteracting system. 



which can be used to parameterize the momentum de- 
pendence. We performed the analytical continuation to 
a;+ = cj + iO by a Pade approximation. Fig. [5] shows 
A(e, uj) together with the integrated intensity A{uj) = 
N~^^^A{e,u)). The analytical continuation is per- 
formed for each spin component, as well as for their av- 
erage. The difference between these results may be re- 
garded as an error of the analytical continuation. 

For n ~ 0.6, i.e., in the Fermi liquid regime [Figs.[2t^a) 
and m^b)], quasiparticle excitations are seen to occur 
around the Fermi energy. Fig. and (d) show the 
spectra for n = 0.8, i.e., for the state with a residual- 
moment. At the higher temperature the dispersion is 
similar to that in the Fermi-liquid regime, but the spec- 
tral features are less sharp than those for n — 0.6. In fact, 
they do not become sharper even for the lowest temper- 
ature of our calculation [Fig. WL^)] ■ We conclude from 
these results, that the phase close to half- filling corre- 
sponds to an incoherent metallic state with a residual 
moment. 

The spectrum shown in Fig. [^d) has very remarkable 
features. Firstly, the broad spectrum crosses the Fermi 



level at the 'momentum' e which is larger than its value 
for the non-interacting case. Secondly, there is an addi- 
tional band above the broad spectrum. We now discuss 
these features in more detail. The effective chemical po- 
tential ^eff = M ~ Rel](0), which is related to the Fermi- 
surface volume if lml](0) = 0, is shown in Fig. [Ijc) to- 
gether with the chemical potentials of the paramagnetic 
and polarized state, /xpM and /.tpM, in the non- interacting 
system [U = J ~ 0). For n < uq, agrees with /ipM, 
indicating that the Luttinger theorem is satisfied. On the 
other hand, it starts to deviate from fipM at n ~ no and 
approaches /xfm- We note that this deviation does not 
imply a violation of the Luttinger theorem, since there is 
no discontinuity in the momentum distribution function 
in this regime. 

The origin of the additional band observed in Fig. [^d) 
can be explained as follows. The weights of the upper and 
lower Hubbard band vary according to n/2. Hence, upon 
hole-doping some weight is transfered from the Hubbard 
bands to energies just above the lower band edge, leading 
to an additional spectral weight S = 1 — n [sH . This addi- 
tional spectrum is due to the hole dynamics. In Fig.[2Id) 
the hole spectrum is separated from the lower Hubbard 
band, while in the coherent Fermi liquid regime the two 
mix. Hence, the appearance of the additional band in 
Fig. [2{d) indicates the incoherence of the holes. A simi- 
lar spectrum has been observed for one hole in the two- 
dimensional t-J model (32j . 

At this point a comment on the possible ground state 
near half-filling is in order. The incoherent metallic state 
has a residual entropy, and therefore an instability is 
expected to take place to lift the degeneracy. Indeed, 
in our calculation with a paramagnetic bath we found 
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FIG. 3: Temperature vs. density phase diagram of the ran- 
dom coupUng t-J model in d = oo for J/t = 0.5. The blue 
dashed curve denotes the boundary of a phase separated re- 
gion. The red dashed curve indicates the transition to a spin- 
glass state. The inset shows the inverse of the spin-glass sus- 
ceptibility Xsg as a function of T. 



a divergence of the charge compressibility dn/d^ below 
T < 0.003, indicating that the paramagnetic solution 
is unstable against phase separation. The temperature 
at which d^/dn changes sign is plotted in Fig. [3] (blue 
dashed curve). Phase separation has also been found in 
the one-dimensional t-J model for large- J [1^. Another 
possibility is magnetic symmetry breaking. Furthermore, 
it is known that at half- filling a spin-glass transition, i.e., 
the breaking of the replica symmetry, occurs [iMl,[il- 
[Tg} , while long-range magnetic order is suppressed due to 
the random distribution of the exchange interaction. To 
estimate the spin- glass transition temperature Tgg, we 
evaluate the spin-glass susceptibility Xag using the ex- 
pression 



Xsg — Ximp/ll J Xiinp)j 



(10) 



derived at half- filling [T^. In the inset of Fig. [3] the tem- 
perature dependence of 1/Xsg is shown. The transition 



temperature Tgg itself is plotted in Fig. |3] (red dashed 
curve). At 71 = 1, we obtain Tsg/ J « 0.147, which 
is consistent with the result of Ref. [ij]- Upon dop- 
ing, Tgg decreases monotonously and reduces to zero at 
n ~ 0.635, a value which is close to no. Therefore the 
incoherent metallic state is actually located in the region 
where the spin-glass phase can be expected to be stable. 
That is, the incoherent metallic state may be stabilized 
if the spin-glass transition is suppressed. Nevertheless, 
the peculiar spectrum in this regime is still physically 
meaningful, since the divergence of Xsg does not affect 
the self-consistency equations. This is the same as in the 
case of the paramagnetic DMFT solution of the Hubbard 
model, which is found in the region where actually the 
antiferromagnetic phase is stable 

In summary, we presented the exact numerical solu- 
tion of the t-J model with random exchange couplings 
in (i = oo. Near half- filling the solution corresponds to 
an incoherent metal with a residual moment. The single- 
particle excitations not only lead to a broad spectrum 
crossing the Fermi level, but to an additional band at 
higher energies. Because of the residual moment, the ef- 
fective chemical potential is located close to the Fermi 
level of the polarized non-interacting system. The ad- 
ditional band is observed only in the non-Fermi liquid 
regime, and is a signature of the incoherence of the holes. 

Finally, we comment on the t-J model with non- 
random couplings in d = 2,3. Eqs. ^ and ([5]) may 
be regarded as a single-site approximation for models in 
Z < oo [131. Indeed, we obtained preliminary results for 
d = 2 and 3 which indicate that the paramagnetic so- 
lution near half-filling is again an incoherent metal with 
a residual moment. The case d < oo, including symme- 
try broken solutions, e.g., the antiferromagnetic or the 
spin-glass state, will be investigated in the future. 
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